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Abstract. Pulsar B0943+10 is well known for its 'B' (burst) mode, characterized by accurately drifting subpulses, 
in contrast to its chaotic 'Q' (quiet) mode. Six new Arecibo observations at 327 MHz with durations of 2+ hours 
each have shed considerable light on the modal dynamics of this pulsar. Of these, three were found to be exclusively 
■ 'B' mode, and three were discovered to exhibit transitions from the 'Q' to the 'B' mode. One of these observations 

QQ ' has permitted us to determine the circulation time of the subbeam carousel in the 'Q' mode for the first time, 

, at some 36.4±0.9 stellar rotation periods. The onset of the 'B' mode is then observed to commence similarly 

fT^ ' in all three observations. The initial circulation time is about 36 periods and relaxes to nearly 38 periods in a 

' roughly exponential fashion with a characteristic time of some 1.2 hours. This is the longest characteristic time 

. ever found in a mode-switching pulsar. Moreover, just after the 'B'-mode onset the pulsar exhibits a symmetrical 

' resolved-double profile form with a somewhat stronger trailing component, but this second component slowly dies 

fH ^ away leaving the usual single 'B'-mode profile with the longitude of the magnetic axis falling at about its trailing 

^ half power point. Thus it would appear that Q-to-B- and B-to-Q-transitions have different characteristic times. 

I ' Some speculations are given on the nature of this slow modal alternation. 
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I. Introduction vestigated. It was shown that the Q-mode onset resulted 

5— i„, , ^ ,,. nir from both abrupt and gradual changes in the intensity 

rrt ' Pulsars are known to lorm very stable mtensity pronles at- r • i- • i i i rn, r i i ■ ^^ 

. . 1 . , 1 1 r • 1- • 1 1 1 All-. 11 or mdividual pulses, ihe presence or slow changes m the 

ter averagmg hundreds or mdividual pulses. Additionally, j t . r ^- i • i r 

, 1 r ,. 11 1 a 1 -1 modal transiormation process having a timescale 01 some 

there are a number or peculiar, so-called mode switch- on • j. j j. j r c ^ i- 

. ,, , 1.1, 11. 1 r., 20 minutes were documented tor the nrst time, 

mg pulsars, which have two stable mtegrated-pronle rr^i • i a -u u ^- ^ j ^ j 

^ ^ ^ , . in c 1 ^"-^^ Arecibo observation was also used to de- 

termine the subbeam configuration responsible for the 
drifting-subpulse pattern for the first time (Deshpande 

,„ , o -r 1 -,no. ri . , , i nnr, ^ Rauklu 1999: 2001, hercaftcr Paper I), and 35-MHz 

(Suleymanova & Izvekova 1984; buleymanova et al 1998, ^ -i . i i ^- ^ ^ ^^ n ^ ^ 

he eafte SIRR") Gauribidanur observations revealed the same 13-mode pat- 

^ , , , , tern at very low frequency (Asgekar & Deshpande 2001, 

In 1992 a nearly 1000-pulse sequence was recorded us- , r. r, tt\ i ui i. r • 

TiiTT 1 • hereafter Paper ii). burthermore, a remarkable set of si- 
mg the Arecibo radio telescope at 430 MHz, and it was .„ j mo tv/ttt i- i ^ 

r . 1 . „„„ , . 1 V. 1 n 1 multaneous 40- and 103-MHz recordings made at the 

found to begin with some 800-1- pulses m the B mode fol- n i i • t-. i- a ; ^i ; . i ■ 

° , „„„ , . . ' 1 rrni r Pushchiuo Rudio Astronomy Observatory suggested in- 
lowed by nearly 200 pulses m its U mode, ihe transfor- , ,,. . 

, 1 11- terestmg correlations between the circulation time, pro- 
mation process from the H to Q mode was thoroughly m- „, „ , . . j i ■ o i c 
1 nle form and emission mode (Rankin, buleymanova & 

Send offprint requests to: joanna.rankinOuvm.edu Deshpande 2003, hereafter Paper HI), but the short dura- 
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shapes. One prominent such pulsar is B0943-I-10. Several 
studies have described the emission characteristics of its 
bright 'B' mode and less intensive (in average) 'Q' mode 
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Fig. 1. Average polarization profile giving the total inten- 
sity (solid line) , total linear (dashed line) and circular po- 
larization (dotted line) for the B-mode of PSR B0943-I-10 
(top panel). Linear polarization angles for the integrated 
profile (solid line) and sufficiently strong individual pulse 
samples (dots) are given as a function of pulse longitude 
(bottom panel). The 327-MHz PS was recorded on 2003 
October 5. 



Table 1. Arecibo 2003 observations at 327 MHz. 



Day 


MJD 


2003 


Mode 


Switch 
(pulse) 


Length 
(pulses) 


1 


52709 


10 March 


Q to B 


2540 


6748 


2 


52711 


12 March 


B 




6809 


3 


52832 


12 July 


Q to B 


5266 


7559 


4 


52840 


20 July 


B 




7274 


5 


52916 


4 October 


Q to B 


1755 


5650 


6 


52917 


5 October 


B 




3024 



It is thus of great interest to further investigate the 
star's modal transformation processes. Observations of 
B0943-I-10 were carried out using the upgraded Arecibo 
radio telescope in 2003 at 327 MHz with this aim. As 
a result, several individual pulse observations have been 
recorded with durations up to 135 minutes. Remarkably, 
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Fig. 2. Average polarization profile and polarization-angle 
distribution as in Fig.[T]for the Q-mode of PSR B0943-I-10. 
This 327-MHz observation was recorded on 2003 October 
4. 

three of the six observations contain PSs of both modes, 
all transitions from the Q to the B mode. One of our main 
questions was whether the transition process between the 
star's two modes is symmetrical. Having heretofore had 
access to only two short records with modal transitions 
(B-Q at 430 MHz in Paper I) and (Q-B at 103/40 MHz in 
Paper III), we now have an excellent opportunity to an- 
swer to this question. §11 discusses our observations, §111 
compares the properties of the two modes, and §IV dis- 
cusses the profile variations incurred in the mode change. 
Fluctuation-spectral variations associated with the modal 
transition are treated in §V, and in §VI we attempt to un- 
derstand where our other B-mode sequences would fall in 
relation to the modal onset. §VII discusses the asymmetry 
between the B-to-Q and Q-to-B transitions processes, and 
§§VIII and IX give brief discussions and summaries of the 
results, respectively. 

II. Observations 

The observations used in our analyses were made using 
the 305-m Arecibo Telescope in Puerto Rico together with 
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Fig. 3. A set of total intensity pulse profiles of B0943+10 
at 327 MHz (2003 March 10), obtained by integration of 
482 consecutive pulses. The first 5 integrations show the 
smooth, broad and weak profiles of the Q-mode. Note the 
unusual profile shape just after the burst-like B-mode on- 
set with the strongest emission coming after the magnetic- 
axis (zero) longitude. The subsequent profiles show the 
gradually decreasing intensity of the second component . 
The broad shape of the total integration in the bottom 
panel shows that Q-like features persist well after the B- 
mode onset. The intensity integrated over the pulse win- 
dow as a function of pulse number is shown in the left-hand 
panel. 



Fig. 4. Another set of profiles of the B0943-I-10 sequence 
in Fig. O obtained here by integration of 27 pulses. The 
abrupt onset of the B mode at about pulse 2540 is asso- 
ciated with a 50% increase in the emission intensity and 
the resumption of subpulse drifting behavior. 



Table 

MHz. 



2. Modal profile properties of B0943-hl0 at 327 



Mode 


Flux 


W 


W 


L/I 


V/I 




(arb.) 


(ms) 


n 


(%) 


(%) 


B 


2? 


20 


7 


30 


15 


Q 


1 


41 


14 


10 


10 



its Gregorian reflector system and new 327-MHz polari- 
metric receiver system. The signals were acquired using 
the Wideband Arecibo Pulsar Processor (WAPP^) dur- 
ing 2003 as detailed in Table 1. The auto- and cross- 
correlation functions of the channel voltages produced 
by receivers connected to orthogonal linearly polarized 
feeds were 3-level or 9-level (days 52916/7) sampled. Upon 
Fourier transforming, 64 or more channels were synthe- 
sized across a 25-MHz bandpass with a 512-/is sampling 
time, providing a resolution of less than a milliperiod. The 
Stokes parameters have been corrected for dispersion, in- 
terstellar Faraday rotation, and various instrumental po- 
larization effects.^ 



III. Characteristics of B0943+10's B and Q modes 

Polarized 327-MHz integrated profiles of B0943-Kl0's B 
and Q modes are given in Figures These observa- 

tions were carried out on 2003 October 5 and 4, respec- 
tively. The top panels give the total intensity (solid), lin- 
ear (dashed) and circular (dotted curve) polarization, and 
the lower panels position-angle (hereafter PA) histograms. 
Table 2 summarizes some of the profile properties. 

The peak intensity of the Q mode is typically 1/3 that 
of its B-mode; however the Q-mode profile is roughly twice 
as wide as its B-mode counterpart. The second component 
of the "pure B"-mode profile is unresolved at 327 MHz — 
similar to the situation at 430 MHz {e.g., Paper I) but 
contrasting with that at 103 MHz, where the second com- 
ponent is 1/3 of the amplitude of main component and well 
resolved, its peak trailing at a distance of 10° (SIRR). 

^ http://www.naic.edu/~wapp 

^ Unfortunately, an error in the WAFF software resulted in 
only one linear polarization being recorded, four-fold redun- 
dant, during the initial 2003 March observations. 
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The lower panels in Figs. show the primary- and 
secondary-polarization modes (hereafter PPM and SPM), 
which form two bands of PA sweep in longitude, sepa- 
rated by 90°. The B mode is dominated by the PPM. The 
PA distribution in Fig. [3 is typical for the Q mode: the 
PPM (lower band) and SPM (upper band) span the entire 
longitude range under the pulse window and the sample 
PAs are almost equally distributed between the two modal 
bands, resulting in a nearly complete depolarization of the 
integrated profile. 

IV. The Q- to B-mode Transitions of B0943+10 
at 327 MHz 

We first discovered the Q- to B-mode transition of 2003 
March 10, and the intensity and profile changes associ- 
ated with this event are shown in Figures O & 01 The cen- 
tral panel of Fig. Ogives a set of 482-pulse averages, five 
before the transition and nine afterward. The transition 
occurred at about pulse 2540 as we will discuss further be- 
low. Note the initial weak, broad, unimodal Q-mode pro- 
file which then gives way to a double B-mode form having 
a more intense trailing component. This unusual trailing 
component then weakens progressively over the next hour 
or so, finally resulting in the usual B-mode profile seen in 
Fig-E — that is, with an asymmetrical single form that has 
longitude zero (that of the magnetic axis) near its trailing 
half-power point. The left-hand panel gives the integrated 
intensity of each average, and it is clear that the B-mode 
onset entails an about 50% increase in the star's inten- 
sity. Fig. 21 gives the same information but in averages 
of only 27 pulses. The breadth of the disorderly Q-mode 
emission is very clear in this figure as is the brightness 
of the occasional intense Q-modc pulses which contribute 
to the averages that are plotted. The two bright features 
between pulses 5500 and 6000 are probably caused by in- 
terference. The lower panels in each figure give the profile 
computed over the full length of the observation and are 
thus identical. Note that this profile is double attesting 
to the strength of the transient second component for a 
substantial period after the modal transition. 

Remarkably, the two other Q- to B-mode transitions 
identified in our observations behave very similarly. The 
one on 12 July occurred near pulse 5266 rather late in 
that day's observation, and a higher level of interference 
on 4 October initially complicated our finding the event 
there at pulse 1755. The left-hand portion of Figure 
shows how the profile changes in all three transitions track 
each other. Here we plot the ratio of the component-height 
maxima A(2/l) as a function of profile number, where 
each profile consists of 256 successive pulses. The three 
transitions are aligned such that profile 1 of each begins 
with its first B-mode pulse. Note how closely the three 
first profiles resemble each other, all with their second 
component peak amplitudes a remarkable 1.75 times their 
first! The diagram shows considerable scatter within the 
first 15 or so profile amplitude ratios, but overall one has 
the impression of an exponential-like behavior. And it is 



Time, hours 
1 2 3 4 5 




-10 10 20 30 40 50 60 
Pulse number (x256) 

Fig. 5. Systematic changes in B0943-|-10s average-profile 
form as a function of time after the B-mode onset. Initially 
the trailing component dominates, but its profile slowly 
evolves to its usual asymptotic form, where the longitude 
of the magnetic axis falls near the trailing half-power point 
and the intensity of the second component decreases to 
1/5 that the first, dominant component. Here, we plot 
the amplitude ratio A(2/l) of the two component peaks, 
separated by 6° longitude, as a function of time in terms 
of the profile number. Each average is comprised of 256 
individual pulses, or 280 seconds (4.7 minutes). Note the 
remarkable similarity of the very first B-mode pulse shape 
with A(2/l) near 1.75 for all three Q- to B-mode transi- 
tion days (hereafter MJD 52709 open triangles, 52832 Xs, 
52916 crosses). Values for the "pure B" days (MJD 52711 
squares, 52840 triangles and 52917 circles) are shifted to 
begin at profile numbers 18, 24 and 22, respectively. They 
clearly show a gradual change of pulse shape and seem to 
represent a later portion of the Q- to B transformation 
process (see also Fig. IIOII . This interpretation is further 
supported by the gradually increasing circulation time of 
all three "pure B" days (see text and Fig.lHJ. The vertical 
dotted line indicates the B-mode onset time. A fit to the 
values (solid curve) has the form 0.17+1.16 exp(— i/r), 
where t is 15.65 256-Pi blocks. 

useful to keep in mind that profile 15 follows the B-mode 
onset by some 70 minutes! 

v. Variation of the 'Carousel' Circulation Time P3 

In Paper HI we noted interesting correlations between the 
forms of the B-mode profiles and the circulation times 
computed from fiuctuation spectra of their constituent 
PSs. The various observations in that paper, though ex- 
tensive, were of inadequate length to see how these varia- 
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Fig. 6. Longitude-resolved fluctuation spectra for the ini- 
tial Q-modc interval of the 2003 March 10 observa- 
tion. The prominent feature in the bottom panel at 
0.0275±0.001 c/Pi corresponds to an interval of 36.4±0.9 
Pi — apparently the first detection of a Q-mode circula- 
tion time. Note also the complete lack of a B-mode "drift" 
feature at about 0.45 c/Pi and the broad, single Q-mode 
profile centered at zero longitude. These MJD 52709 spec- 
tra were computed using the first 2048 pulses, whereas the 
B-mode recommences at about pulse 2540. 

tions came about. Moreover, it was unclear whether sub- 
beam circulation persisted in the Q mode; we found some 
hints in Paper I that it did do, but no Q-mode PS then 
in our possession exhibited a fluctuation feature demon- 
strating such circulation. 

We were thus surprised and gratified when the Q-mode 
PS preceding the B-mode onset of 2003 March 10 ex- 
hibited a strong low frequency modulation feature. The 
longitude-resolved fluctuation (hereafter Irf) spectra are 
given in Figure El They are computed from the first 2048 
pulses in the PS having a B-mode onset at pulse 2540. A 
"boxcar" smoothing of 3 samples has been applied to the 
PS to enhance the S/N of the spectra. The Q-mode profile 
is clearly recognizable in the left-hand panel and the inte- 
gral spectrum is shown in the lower panel where a strong 
feature at some 0.025 cycles/period (hereafter c/Pi) is evi- 
dent. Careful measurement of the frequency of this feature 
resuhs in a value of 0.0275±0.001 c/Pi. This is the first 
time that a fluctuation feature corresponding to the cir- 
culation time has been seen in our meter-wavelength ob- 
servations, though such a feature has at times been seen 
in the decameter PSs of Paper II. The circulation time 
associated with this fluctuation frequency is then some 
36.4±0.9 Pi — a somewhat lower value that usually found 
for B-mode circulation times. 



Average profile Fluctuation spectra 

2x10'' (longitude-resolved) 




0.1 0.2 0.3 0.4 
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Fig. 7. Longitude-resolved fluctuation spectra computed 
from 256 pulses just after the B-mode onset of the 
2003 March 10 observation. Here the "drift" feature fre- 
quency is 0.452 c/Pi, which corresponds to a 20-subbeam 
"carousel" circulation time of 36.47±0.03 Pi — ^just the 
same as that determined for the Q-modc in Fig. El 

Subbeam circulation times can more readily be com- 
puted for B-mode PSs from the strong roughly 0.46 c/Pi 
subpulse-modulation feature always associated with this 
behavior. We further know that this response is a first- 
order alias of the true fluctuation frequency of about 0.54 
c/Pi — a matter which has been demonstrated extensively 
using the harmonic-resolved fluctuation spectrum {e.g., 
Paper I). Finally, every analysis of the star's B-mode prop- 
erties has resulted in the conclusion that the observed sub- 
pulse drift is produced by a circulating pattern of just 20 
subbeams. Thus, the B-mode circulation time can be com- 
puted as P3=20P3(true)=20/[l - /3(obs)]. 

We give an example of such a computation in Figured 
for the 256 pulses immediately after the B-mode onset 
of the 2003 March 10 observation at pulse 2540. Here 
the highly unusual profile just following the B-mode on- 
set is given in the left-hand panel, and the bottom panel 
shows the integral spectrum associated with the drifting- 
subpulse modulation. The feature frequency is here some 
0.452±0.003 c/Pi — one of the smallest such values ever 
measured for the star's B mode. 

Generally, it was possible to measure a circulation time 
for each 256-pulse interval following the respective B-mode 
onsets in each of the three observations. Not only do the 
P3 values of these three transitions behave similarly as 
a function of time, but they too suggest an exponential 
dependence of the form [1 — cxp(— i/r)]. The behavior of 
all three observations is plotted in the left-hand part of 
Figure (HI where the three observations are depicted with 
triangle, -I- and x symbols. Note as well that several other 
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Fig. 8. Subbeam "carousel" circulation time P3 as a func- 
tion of time from the B-mode onset at intervals of 256 Pi; 
1 hour is equal to 12.8 such units. Circulation times for the 
three Q- to B-mode transition observations are aligned at 
the origin (vertical dotted line) and fitted with the expo- 
nential curve (solid line). The other three "pure B"-mode 
observations could be fitted to this curve by shifting along 
the time axis at 18, 24 and 22 (256- Pi units), respectively. 
The characteristic time is some 4000 pulses or 1.22 hours. 
It is also about 107 circulation times. All three days with 
"pure B" -profiles show a gradual P3 increase from some 
37.3 to about 37.8 Pi during the observational session and 
thus seem to be a later part of the Q- to-B transformation 
process. Four of the Irf spectra exhibited double peaks and 
P3S calculated for these weaker precursor peaks are sep- 
arately indicated (open circles) but not included in the 
fitting procedure (see text). 



of the circulation time over their duration. Of course, we 
now understand that a Q- to B-mode transition had to oc- 
cur somewhat before the beginning of our data-acquisition 
in each of these cases — the only question is just how long 
before. 

The profile and circulation-time variations give us an 
indication of just where the "pure B"-mode observations 
fall in relation to the B-mode onset. If we proceed on 
the strong appearance that these quantities do exhibit a 
[1 — exp(— t/r)] functional dependence, then we can lo- 
cate the various average profiles and their corresponding 
P3 values appropriately and then partially test whether 
the procedure was correct. Indeed, we find again that the 
two measures vary together such that the three days can 
be interpolated at a time corresponding to the 18th, 24th 
and 22nd 256-pulse average after the B-mode onset, re- 
spectively. We have plotted the A(2/l) values in Fig. |51 
and the P3 measurements in Fig. |S1 as is also described in 
their captions. 

While the above interpolation procedure must be re- 
garded as approximate, it also cannot be wildly wrong. 
The above circulation-time values following the B-mode 
onset were fitted to a [1 — exp(— t/r)] function, which is 
shown by the continuous curve in Fig. |S| The character- 
istic time r is then about 73 minutes or 4000 Pi. It is 
also some 107 circulation times. In terms of the light-time 
scale of the pulsar's magnetosphere, to say nothing of the 
acceleration regions within its polar flux tube, this is a 
very long time indeed. It is hard to imagine what physical 
processes associated with pulsar emission could have such 
a long time scale. 

Finally, in Figure El we plot the component peak- 
amplitude ratio A(l/2) versus the subbeam-circulation 
time P3 for all of six observations as well as some fur- 
ther observations from Paper III (open square symbols). 
The values are well fitted by an inverse relation. Overall, 
it is clear that the subbeam-circulation time and profile 
form are closely related, such that one could reliably be 
estimated from the other. 



of the initial circulation-time values fall at just more than 
36 Pi , thus agreeing well with the more poorly determined 
Q-mode value from Fig. |^ 

Remarkably, we can then apparently conclude that the 
initial B-mode, subpulse-drift associated circulation time 
represents a continuation of that associated with the very 
disorderly emission of the star's Q mode. 

Vl. Locating the "Pure B"-mode Observations 
Within the Transition Recovery 

We then turned our attention to the analysis of the three 
long PSs in which only B-mode emission was observed, 
those recorded on 12 March, 20 July and 5 October, re- 
spectively. Each PS exhibited two properties which sur- 
prised us initially: a) there was a perceptible tendency 
for the amplitude ratio A(2/l) to decrease, and b) we also 
found strong evidence in two of the three for a lengthening 



Vll. Asymmetry in the Mode-Changing Process 

We see good evidence from the foregoing analyses that 
the modal-transformation process of the emitting regime 
in pulsar B0943-I-10 is not symmetrical in time. Thus 
far we have only managed to observe one B-to-Q-modal 
transition — this the remarkable 18-min (986- Pi) Arecibo 
430-MHz observation obtained in 1992, which was thor- 
oughly investigated in both SIRR and Paper I. 

Comparison of the evolution of the pulse shape and 
subbeam-carousel rotation time P3 for the respective Q-to- 
B (2003) and B-to-Q (1992) transitions reveals significant 
differences — i.e., in the latter case these processes occur 
much more quickly. 

To investigate the pulse-shape evolution we divided 
the pulse window onto two subwindows Wl and W2 at 
negative (~8°, 0) and positive longitudes (0, +8°). As was 
shown earlier at 430 MHz, the maxima of the B-mode and 



< 4 
o 



B 3 



CO 

to 



o 



o 

O 1 



Rankin & Suleymanova: Pulsar 0943+10's 'Drifting'-Subpulse Emission IV. 
I I I I I I I I I I I I I H Time, hours 

. - -2-10123 



□ 



□ 



+ 

X 



+ A 



A A + 



36 36.5 37 37.5 

Subbeam circulation time, P1 



38 



Q-B transition 
327 MHz (2003) 



5 3 — 



g 
13 



S 2 — 
c 

B 




B-Q transition 
430 IVIHz (1992) 



I I I I I I 
150 200 



250 



Hundreds of rotation periods 



Fig. 9. Component amplitude ratio A(l/2) as a function 
of subbeam circulation time P3 at 327 MHz. The 103-MHz 
values (open square symbols) were taken from Table 1 of 
Paper III. 

Q-mode profiles are separated in longitude at a distance 
nearly equal to the half-power width of the B-mode pro- 
file. This means that a longitude on the trailing edge of 
the B-mode profile corresponding to the half-power point 
can be taken as a zero longitude, equally remote from 
regions where the dominance of each mode is strongest. 
Temporal variations of the intensity ratio in these two 
subwindows W1/W2 are presented in Figure^] The in- 
tensity of the subwindow for each individual pulse is ob- 
tained by summing the amplitudes of all the bins for a 
given subwindow — besides additionally averaging over 100 
adjacent pulses to diminish the influence of the noise and 
possible scintillation fluctuations. Note that the W1/W2 
ratio changes from 2.8 to 0.8 during the 18-min B-to-Q 
profile-shape transformation; whereas in the opposite Q- 
to-B transition the same change in pulse shape takes about 
2.5 hours! Similarly, the P3 value at the B-mode recom- 
mencement is just over 36 Pi, while at the Q-mode onset 
it was about 37.35 Pi. One can suggest that a sharp de- 
crease of the circulation time down to 36 Pi occurs after Q- 
mode commencement, without further significant changes 
of the pulse shape. Note that the Wl /W2 ratio is nearly 
constant (W1/W2 is randomly distributed in the range 
0.6-1.0) during the initial 46-min (2540 pulses) Q-mode 
interval of 10 March and the 96-min (5266 pulses) Q-mode 
PS of 12 July 

VIII. Discussion 

Pulsar B0943-|-10's circulation time P3 in Fig. |H1 gradu- 
ally changes in the range from about 36 Pi to at least 



Fig. 10. Intensity ratio of two emission sub-windows for 
negative longitudes Wl and positive longitudes W2 as a 
function of time following the B-mode onset in units of 
the group number of successive 100-pulse averages. The 
pulse shape is apparently stable during the 1.5 hours of 
Q-mode emission, whereas it continuously changes after 
the B-mode recommencement. The mode transition times 
are shown by vertical dashed lines. The shift of the 1992 
B- to Q-mode observation (430 MHz) along time axis is ar- 
bitrary. This figure clearly shows the different pulse-shape 
transformation rate for Q-to-B- (slow) and B-to-Q-mode 
(fast) transitions in B0943-I-10. 

37.8 Pi following the onset of its B mode. The difference 
(1.8/36=0.05) may be as large as 5%. 

Drift-rate variations are well known to occur in a num- 
ber of pulsars, and we now see that B0943+10 is an- 
other member of this group. Some stars such as B0031-07 
{e.g., Wright & Fowler 1981a; Vivekanand & Joshi 1997), 
B1944+17 {e.g., Deich et al 1986), B2303-f'30 (Redman 
et al 2005) and B2319-F60 {e.g., Wright & Fowler 1981b) 
exhibit several discrete drift-rate "states" that are usu- 
ally designated as modes. Others show what appear to be 
continuous drift-rate changes from rarely {i.e., B0809-I-74; 
van Leeuwen et al 2002) to very frequently (B0826~34 
and B2016-I-28)— and in B0031-07 and B2303-f-30 a regu- 
lar cycle has been detected following occasional "events" 
wherein "null" pulses are often observed to accompany 
drift-rate variations. 

While the discrete drift rates can change by several 
factors of two, the amplitude of observed continuous drift- 
rate variations is less than 8% as in B0943-I-10, and 
their characteristic times some several minutes — i.e., sig- 
nificantly less than we observe in B0943+10. Moreover, 
B0943-f 10 surely provides the best example of drift-rate 
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changes which correlate with large changes in the profile 
form. In most other such pulsars drift-rate variations are 
associated only with smaller and less orderly changes in 
the amplitude, shape and/or phase of the pulse profile. 
Systematic variations of the drift rate in B0826-34 (and 
probably B2303+30) apparently cross alias boundaries 
and thus produce changes in the drift direction (Esamdin 
et al 2004; Gupta, Gil, Kijak, & Sendyk 2004). An even 
stronger similarity with B0943+10 is the discovery by 
Esamdin et al that B0826-34's long intervals of "null" 
pulses are in fact a weak Q-like mode, which has a very 
different profile from that of the strong, drifting mode (see 
their fig. 2). This weak mode lasts longer than the stronger 
mode, and its profile is also found to be stable in intensity 
and shape. Moreover, both stars are now known to have 
closely aligned geometries. 

In B0943+10 it is notable that null pulses have not 
been detected despite its prominent drift-rate variations. 
In its bright B mode an upper limit on the number of 
nulls might be placed at around 0.1%, but in the chaotic 
Q mode the intensity range of individual pulses is so great 
that no sensible discrimination can yet be made between 
possible nulls and weak pulses. 

An effort to understand drift-rate variations in the 
context of revisions to the polar-cap accelerator theory 
of Ruderman & Sutherland (1975) has recently been pub- 
lished by Gil, Melichidze & Geppert (2003) . These authors 
attribute such variations to the shielding effect of iron ions 
whose flow is sensitive to the surface temperature of the 
polar cap, which in turn is heated by the bombardment 
of ultra-relativistic charged particles originating from pair 
production. The model implies that the cooling process of 
the star's surface below the mean temperature leads to a 
faster than mean drift rate, while the reverse process of 
heating results in slower than mean drift rates. The drift- 
rate and direction changes of B0826-34 at 325 MHz were 
discussed by Gupta, Gil et al (2004). A typical time scale 
for one cycle of drift-rate reversal is about 100 stellar rota- 
tion periods, or nearly 3 min for this pulsar. In B0826 34 
drift-rate increases (corresponding to cooling below the 
mean temperatiue) are found to occur more slowly and 
smoothly than drift-rate decreases. 

In B0943-I-10 the cycle is just the opposite: the slowly 
decreasing B-mode drift rate (corresponding ostensibly to 
heating above the mean temperature) is the slower pro- 
cess. Thus, this mechanism would appear to be backward 
for B0943-hl0, contrary to B0826-34. Nevertheless, ft is 
clear that the circulation-time increases and decreases of 
the subbeam carousel have different characteristic times 
for the two pulsars. The current theory predicts very 
short cooling time scales (a few /zsec for magnetic fields 
strengths of a few 10^^ G; Gil et al 2003) that could explain 
microstructure of individual pulses or short nulls. That 
any physical mechanism could cause multi-hour-scale vari- 
ations in the polar cap temperature in pulsars is still very 
much open to question. 

As mentioned above, a number of pulsars exhibit 
drift-rate variations in a range up to about 8% like 



B0943-I-10. None of the others, however, demonstrate 
B0943-|-10's pronounced and correlated changes in profile 
form. Though we cannot be sure what agent or mecha- 
nism is causing these profile variations, we are attracted 
by the analogy to the average pulse-shape changes in two 
well studied "mode switchers", B0329+54 and B1237+25. 
Core emission is prominent in these two stars, surely in 
substantial part because our sightline passes much closer 
to their respective magnetic axes, but in both cases clear 
modal variations are also seen in the characteristics of 
their conal components. A tendency for profile modes 
to represent more and less symmetrical configurations of 
emission (about the longitude of the magnetic axis) has 
long been noted {e.g., Rankin 1986), and the contrast- 
ing properties of the modes in these two stars is partic- 
ularly striking in this respect: both exhibit roughly sym- 
metrical profiles in their respective "normal" modes and 
then strong leading/ trailing emission asymmetries accom- 
panied by enhanced core activity in their "abnormal" 
modes. In B0329-I-54 at meter wavelengths {e.g., some 
111-606 MHz), mode changes are marked by alterations in 
the positions and intensities of the core and leading conal 
components (Gangadhara & Gupta 2001; Suleymanova & 
Pugachov 2002), and a very similar pattern is observed in 
B1237-^25 (Bartel et al 1982; Srostlik & Rankin 2005) as 
well apparently as in B0355-I-54 in the course of its slow 
modal variations (Morris et al 1980). 

In B0943-I-10 [as in other conal single (S^) stars] our 
sightline cuts across the periphery of its conal emission 
pattern, so we can have no direct indicator of its core ra- 
diation intensity. Nevertheless, we can hypothesize that 
B0943-f 10 belongs to the group of pulsars with variable 
core activity. In B0329+54 and B1237-I-25, core-active 
intervals are observed as sudden eruptions, whereas in 
B0943-I-10, this activity may increase gradually over sev- 
eral hours. As a result the star's symmetrical profile (seen 
in its Q mode and also just after B-mode onset) evolves 
to its usual, highly asymmetrical form several hours into 
each B-mode interval. It is also possible that the pro- 
file change at B-mode recommencement from "single" to 
"double resolved" entails a conal emission-height increase 
that is prompted by core activity (indeed, a possibility 
worth exploring for B0826-34 as well). 

Recent XMM-Newton observations by Zhang et al 
(2005) have detected X-ray emission from B0943-I-10, and 
several current models suggest that such radiation is basi- 
cally thermal, emitted from a polar cap surface heated by 
back-flowing plasma produced by the same pair-processes 
that are responsible for the star's radio emission. These 
authors conclude that "within the thermal interpretation, 
the X-ray radiation is emitted from a heated area much 
smaller than the conventional polar cap area" . This result, 
however, was based on a 20-ks observation — one surely 
long enough to involve one or another (or a mixture) of B 
and Q-mode intervals. Increased core activity (associated 
with one of the modes) might alter the polar-cap surface 
temperature and distribution (if the cooling time-scales 
are indeed long enough), affecting ExB and/or the conal 
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beam symmetry as mentioned above. The above result 
takes no account of this and is thus of uncertain signifi- 
cance. 

If, however, the "feedback" model of Wright (2003) 
is substantially correct, the situation may be much more 
complex. Hero the "carousel" circulation time would de- 
pend on the "connectivity" between the polar cap and 
"outer gap" accelerators, and only indirectly on the ExB 
drift. In this aurora-like scenario, the variation in both ro- 
tation speed and profile asymmetry would be produced by 
the inevitably asymmetric interaction between the mag- 
netic poles of an inclined rotating star. Even in a perfect 
dipolar geometry, the slowly circulating emission columns 
at either pole would not remain on matching fieldlines 
and their communication paths would gradually become 
"twisted" in a manner not unlike the earlier suggestion of 
Paper III. 

IX. Summary and conclusions 

B0943+10's subpulse-drift rate is observed to decrease ex- 
ponentially by some 5% during the 4+ hours after B-modc 
onset with a characteristic time of some 73 min — and a 
very similar behavior is seen in all three 327-MIIz ob- 
servations. These drift-rate changes are accompanied by 
changes in the integrated pulse shape. Remarkably, the 
trailing component brightens at B-modc onset, making 
the overall profile a conal double, and for the first 20 
mins or so at 327 MHz it outshines the usual leading 
emission feature. Overall, the amplitude ratio of the two 
B-mode profile components varies significantly from 1.75 
to 0.2 during the same 4+ hours. Apparently this slowly 
decreasing drift rate eventually falls below some critical 
value, and is replaced by the disorderly drift pattern of 
the alternative Q emission mode. This Q mode is charac- 
terized by a broad Gaussian-shaped profile that is stable 
over several hours, a complete absence of regular subpulse 
drift, and a diminuation to about half in emission inten- 
sity. This highly organized process of gradually decreas- 
ing subpulse drift rate accompanied by systematic profile 
variations has no simple parallel in other radio pulsars. 
In particular, no other known pulsar exhibits systematic 
variations on a time scale even approaching the 1.2 hours 
seen in B0943+10. 

Our main results can be summarized as follows: 

— Six observations at 327 MHz with a total duration of 
more than 11 hours were carried out in 2003 in order 
to search for modal transitions. 

— Three of the days exhibited Q- to B-modc transitions 
and three others were found to have long B-mode se- 
quences. 

— All three Q- to B-modc transitions appear to behave 
very similarly, qualitatively and quantitatively. 

— A Q-mode circulation-time feature was detected for 
the first time, and its value is some 36.4±0.9 Pi. 

— The B-mode circulation time is found to be some 36.2 
Pi just after its onset. 



— The B-mode circulation time exhibits an exponential 
relaxation to an asymptotic value of perhaps 37.8 Pi 
with a characteristic time of 1.2 hours. 

— The early B-mode profile is found to be double — with a 
somewhat larger second component — which then grad- 
ually decreases over a several hour period, such that 
the magnetic-axis longitude falls at the trailing half- 
power point of the asymptotic B-mode profile. 

— The symmetrical to asymmetrical subbcam-carouscl 
visibility — demonstrated by the dramatic B-mode pro- 
file variations — points to a complex physical origin re- 
lated to the "absorption" phenomenon. 

— Q-to-B-mode changes in pulsar B0943-I-10 again have a 
rapid and very slow aspect. The 1.2-hour characteristic 
time may be the longest orderly behavior known to be 
associated with pulsar emission. 
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